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INTRODUCTION 

The development of two-phase flow research under 
reduced and microgravity conditions is prompted by a 
wide range of space applications, such as, thermal en- 
ergy and power generation, propulsion, cryogenic stor- 
age and long duration life support systems; necessary 
for programs, such as NASA’s Human Exploration for 
the Development of Space (HEDS). Study of gas- 
liquid flows in reduced gravity is important for the 
design of two-phase thermal control systems, intended 
to replace conventional pumped liquid loops 1 . The 
main advantage of a two-phase thermal control system 
is its reduced weight, which is of utmost importance 
since any weight saving for payload transfer to orbit 


signifies lower launch costs. 

Under normal gravity conditions, when gas is in- 
jected through an orifice into a quiescent liquid me- 
dium, the bubble grows and detaches quite readily due 
to the buoyancy force. Under reduced gravity condi- 
tions, the detaching role of the buoyancy force is sig- 
nificantly diminished, giving rise to uncontrollably 
larger bubbles than those obtained in normal gravity. 
Consequently, another bubble detaching force is re- 
quired in order to control bubble size and frequency of 
formation. A practical solution is to use flowin^liquid 
and utilize its drag force for bubble detachment ' . 

Two configurations generally considered for bub- 
ble dispersion in a flowing liquid are the co-flow and 
the cross-flow geometry. In the co-flow configuration, 
the dispersed phase is introduced through a nozzle in 
the same direction with the liquid flow; whereas in the 
cross-flow geometry, gas is injected perpendicular to 
the direction of liquid flow. Of the three major two- 
phase flow patterns, namely bubble, slug and annular , 
only bubble and annular flows are used in space based 
systems 6 . In this work, we are concentrating on the 
bubbly flow regime. 

For normal ’ 8 and reduced gravity 9 ' °, bubble gen- 
eration in a quiescent liquid has been extensively 
studied. Reduced gravity bubble formation in the 
cross-flow configuration has been recently reported by 
several investigators 4, lI ’ n . On the other hand, the co- 
flow configuration has only been considered in normal 
gravity in order to observe the effect of liquid velocity 
on bubble detachment 12, l3 . Under reduced gravity 
conditions bubble generation and resulting two-phase 
flow by multiple nozzle injection along the periphery 
of the flow conduit has been reported by several in- 
vestigators 5 ’ 6 . With multiple nozzle injection, due to 
unpredictable coalescence of adjacent bubbles, it is 
difficult to control accurately the void fraction of the 


ensuing two-phase flow. A better alternative is con- 
trolled bubble generation via single nozzle injection. 

In this work, we investigate bubble generation by 
gas injection via a single nozzle in a co- and cross- 
flow system. Experiments using air and water are per- 
formed in parabolic flight aboard the modified DC-9 
Reduced Gravity Research Aircraft at NASA Lewis 
Research Center. Effects of surrounding liquid velocity 
and two-phase flow conduit geometry on the bubble 
diameter and the associated void fraction are investi- 
gated. For the co-flow geometry, we have also devel- 
oped a theoretical model, based on an overall balance 
of forces acting on the bubble. Predictions of bubble 
diameter and formation frequency using the present 
model show good agreement with our experimental 
results. 

EXPERIMENTAL STUDY 

In reduced gravity, the bubble diameter (D B ) is 
dependent on fluid properties, flow geometry and flow 
conditions. 

D B = f ( Qd, Qc, Dp, Dn, a, li c , fid, Pc, Pd ) 
where Q d is the volumetric gas flow rate, Q c is the 
volumetric liquid flow rate, D P is the pipe diameter, D n 
is the nozzle diameter, a is the surface tension, |i c and 
p c arc the dynamic viscosity and density of the liquid 
phase, while ^ and p d arc the dynamic viscosity and 
density of the gas phase. In this study, we are investi- 
gating the effect of flow conditions (Q d , Q c ) and flow 
geometry (D P , D\ ) on bubble diameter (Dg). 

Apparatus and Methods 

Our experiments are conducted aboard the modi- 
fied DC-9 Reduced Gravity Research Aircraft which 
provides the investigator with 20 seconds of 0.01 g 
reduced gravity environment. Out of this time period 
an estimated 15 seconds is allotted for acquisition of 
experimental data. The reproducibility of data between 
two consecutive trajectories is within ±5%. The pres- 
ent two-phase flow experiments are performed with an 
air-water system using three different sets of pipe di- 
ameters (D p = 1.27 cm, 1.9 cm and 2.54 cm). In addi- 
tion, two different ratios of nozzle to pipe diameters 
are considered (D n * = = 0.1 and 0.2). Depending 

on the two-phase flow pipe diameter, superficial gas 
(U G s = 4Q d /7tD p 2 ) and liquid velocities (Uls = 
4Q c /7iD p 2 ) are varied from 8 to 70 cm/s. In general, for 
each data point acquired, two dive trajectories are exe- 
cuted. 

The co- flow test section (figure 1) consists of a 
Plexiglas pipe, which acts as the two-phase flow con- 
duit. 
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Fig 1 : Co-flow experimental test section 


A tee branch fitting is mounted on the inlet side of the 
pipe. Dry and Altered air is injected through a stainless 
steel tube, which acts as the gas injection nozzle and 
protrudes into the pipe. Distilled water is introduced 
through the remaining port of the tee branch fitting. 
The water and air mixing region is surrounded by a 
Plexiglas visual rectangular box filled with water, 
which eliminates optical distortion of the generated 
bubble. 

The cross flow test section, shown in figure 2, is 
machined as a tee section from a rectangular piece of 
Plexiglas stock. Two orthogonally positioned, equal 
diameter holes are bored into the Plexiglas tee, one 
merging into the other. Through one of these holes, air 
is injected via a stainless steel tube. The other hole acts 
as the water inlet tube as well as the two-phase flow 
conduit. 



Fig 2: Cross-flow experimental test section 


The complete experimental test section assembly is 
integrated within the Lcarjet Two-phase Flow Appa- 
ratus, developed and described by McQuillen and 
Neumann 14 . The bubble generation process is recorded 
with a high speed video camcorder (250 frames/sec). 
For every experimental run, a new batch of distilled 
water is used to minimize contamination at the bubble 
surface. 


Experimental bubble diameter is obtained from the 
flight experiment video by using THIN 2.0® and 
OPTIMAS 5.1 c image acquisition and processing 
software packages. The geometrically averaged bubble 
diameter for each of three consecutively detached bub- 
bles in the vicinity of the gas injection nozzle is first 
calculated. The bubble diameter reported in this work 
is the arithmetic average of these three values. The 
standard deviation for bubble diameter measurement is 
within ± 2.5% of the mean diameter value. Uncertainty 
errors in flow velocity measurement have an upper 
limit of ± 5% of the obtained value. A more elaborate 
discussion on experimental procedure and data acqui- 
sition is given by Pais 15 . 


Co-flow Configuration Results 

The important role played by the flowing liquid on 
bubble detachment in reduced gravity is shown in fig- 
ure 3. Experimental data presented in this figure is 
obtained for a fixed pipe diameter of 1.9 cm and a gas 
flow rale of 51 cm /s (cc/s). Two sets of nozzle di- 
ameters (D n = 0.19 cm and 0.38 cm) are used. For both 
nozzle diameters, as superficial liquid velocity is in- 
creased, the bubble diameter decreases. The drag in- 
duced by the surrounding liquid flow results in faster 
detachment of the bubble, thereby decreasing its size. 
It is further observed that bubble diameter increases 
with nozzle diameter, which can be explained from the 
fact that with a larger injection nozzle more gas is fed 


into the bubble. 
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Fig 3: Effect of nozzle diameter and superficial liquid 
velocity on bubble diameter. Fixed Qd = 51 cc/s and 
Dp = 1.9 cm. 


Variation of bubble size with increasing volumetric 
gas flow rale and pipe diameter is displayed in figure 
4. The volumetric liquid flow rate and the nozzle di- 
ameter are kept constant at Q c = 40 cc/s and D>, = 0.19 
cm. Experiments are conducted with three different 
pipe diameters, D P = 1.27, 1.9 and 2.54 cm by varying 
the gas flow rate from 16 to 40 cc/s. It is obvious that 
bubble diameter increases with increasing gas flux and 
pipe diameter. At a fixed Q c , an increase in pipe di- 
ameter implies a reduction in co-flowing liquid veloc- 
ity and therefore liquid drag. Larger bubbles are 
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formed when drag induced by the flowing liquid is 


reduced. 
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Fig 4: Variation of bubble diameter with volumetric 
gas flow rate and flow conduit diameter. Fixed Qc = 
40 cc/s and D N = 0. 1 9 cm. 


Figure 5 displays variation of void fraction (V f ) 
with volumetric gas and liquid flow rates. Void frac- 
tion is defined as the ratio of volume occupied by the 
gas phase to total volume of fluid within a given sec- 
tion of the two phase flow conduit. Mathematically, 
for a single bubble in the bubbly flow regime, this re- 
lationship can be written as Vf = 2D B 3 /3D P 2 A, where A 
is the distance between the front of the detached bub- 
ble and the front of the previously detached bubble. In 
the bubbly flow regime the maximum value of the void 
fraction is 2/3 (D B = D P ), beyond which formation of 
Taylor bubbles ( slugs ) occurs. In the slug flow re- 
gime, void fraction can exceed the value of 2/3, its 
upper limit being 1 . 
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Fig 5: Dependence of void fraction on gas and liquid 
flow rate with respect to variation in nozzle diameter. 
D P is fixed at 1.9 cm. 

Void fraction, which is directly proportional to the 
bubble diameter, increases with gas flow rate and noz- 
zle diameter, while decreasing with liquid flow rate. 
For the data presented in this graph, the pipe diameter 
is fixed at 1.9 cm. It is observed that at D N * = 0.1 (D N 
= 0.19 cm), for two different liquid flow rates, Q c = 35 
and 68 cc/s, the void fraction is less than 2/3, indicat- 
ing occurrence of the bubbly flow regime. On the other 
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hand, for Q c = 35 cc/s and D N * = 0.2 (D N = 0.38 cm), 
the void fraction exceeds the value of 2/3 for Q d = 29 
and 33 cc/s, suggesting formation of Taylor bubbles at 
such flow conditions. 


Cross Flow Configuration Results 

In the cross-flow configuration, bubble diameter as 
a function of superficial liquid velocity and nozzle 
diameter is shown in figure 6. This plot displays data 
taken with the 1.27 diameter test section at a constant 
gas flow rate of 44 cc/s. For two different nozzle di- 
ameters (D n = 0.127 and 0.254 cm), the superficial 
liquid velocity (U L s) is varied from 20 to 60 cm/s. It is 
observed that the bubble diameter decreases with in- 
creasing superficial liquid velocity for a given gas in- 
jection geometry at a constant gas flow rate. Further- 
more, it is noted that the bubble diameter increases 
with increasing nozzle diameter. The trends are similar 
to those displayed by the co-flow configuration. 



Fig 6: Effect of superficial liquid velocity and nozzle 
diameter on bubble size for a fixed Q d = 44 cm /s and 
D P = 1.27 cm. 
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Fig 7: Variation of bubble diameter with respect to gas 
flow rate and pipe diameter at a fixed Q c = 40 cm 3 /s 
and Dn = 0.19 cm. 


Figure 7 shows the variation of bubble diameter 
with volumetric gas flow rate and pipe diameter at a 
constant liquid flow rate of 40 cc/s. For acquiring this 
data, three distinct volumetric gas flow rates are used, 
namely 16, 24 and 40 cc/s. Experiments are performed 
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using the 1.27, 1.9 and 2.54 cm diameter test sections, 
keeping the nozzle diameter constant at 0.19 cm. Note 
that bubble size grows with increasing gas flow rate 
and pipe diameter at a constant liquid flow rate. Figure 
4 and Figure 7 display bubble diameters at identical 
Bow geometry and conditions for co and cross flow 
configurations. It is interesting to observe that at simi- 
lar flow conditions and geometry, somewhat larger 
bubbles are generated by using the cross-flow configu- 
ration rather than the co-flow configuration. 

A plot of void fraction (V f ) values as a function of 
volumetric gas and liquid flow rates is presented in 
figure 8. This experimental data is obtained using the 
1.9 cm diameter test section for two different nozzle 
diameters, D N = 0.19 and 0.38 cm. For two different 
liquid flow conditions (Q c = 51 and 68 cc/s), the gas 
flow rate is varied from 16 to 49 cc/s. In the cross-flow 
configuration, analogous to the co-flow geometry, the 
void fraction increases with increasing gas flow rate 
and decreases with increasing surrounding liquid flow. 
Furthermore, it is observed that V f increases with gas 
injection nozzle diameter. 



volumetric gas flow rale, Qd (cm 3 /s) 


Fig. 8: Effect of volumetric gas and liquid flow rates 
on void fraction with respect to variation in nozzle 
diameter for constant pipe diameter of 1 .9 cm. 

Comparison of Co- and Cross-flow Configuration 

Bubble diameters and void fraction obtained using 
the co-flow system are compared with corresponding 
values obtained in the cross-flow geometry. The data 
displayed in figure 9, is obtained using a 1.9 cm di- 
ameter test section at constant liquid flow rate of 68 
cc/s with a 0.38 cm nozzle diameter. The volumetric 
gas flow rate is varied from 21 to 70 cc/s. It is ob- 
served that at similar values of gas and liquid flow 
rates as well as similar nozzle and pipe diameters, 
bubbles generated by using the cross-flow configura- 
tion are slightly larger in size relative to those obtained 
in co-flow geometry. Therefore, the void fraction of 
the resulting two-phase flow also follows a similar 
trend. 
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Fig 9: Comparison of bubble diameter and void frac- 
tion for co and cross-flow configuration. Fixed Q c = 68 
cc/s, Dp = 1.9 cm and D N = 0.38 cm. 

THEORETICAL STUDY 

In parallel to the experimental work, we have de- 
veloped a theoretical model to describe the bubble 
detachment process from which wc obtain the de- 
tached bubble diameter and formation time in a co- 
flow configuration, as shown in figure 10. 





Fig. 10: Schematic of co-flow configuration. 

The current model, which is valid for both normal 
and reduced gravity conditions at constant gas flux, is 
employed to investigate single bubble generation in the 
dynamic (Q d = 1 - 1000 cc/s) and bubbly flow regime 
(D b < Dp). The focus of the model is to identify the 
important forces involved in the process of bubble 
formation and their role on bubble detachment. The 
bubble shape is assumed spherical throughout the for- 
mation process. Hence, for constant flow conditions, 
the rate of change of bubble volume (V B ) is given as d 
V B /d t = Q d = constant. Various important forces in- 
volved in the bubble generation process are: 

F B - V B (p c -p d )s 

F„ = JtD N o 

F =p Q ^- 

M Pd (rc/4)D * 

Fd — S D C DW — p c U c ff A cfr 
F, =^[p<V B f-+PcC MC V B U c(T 
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where F b , F 0 , F M , F D , F! are respectively the buoyancy, 
surface tension, momentum flux, drag and inertia 
forces acting on the bubble. Cow is drag coefficient of 
the bubble in presence of confining pipe walls; ds/dt is 
velocity of the bubble center away from the nozzle tip; 
Ueff is the relative velocity of bubble center w. r. t. 
superficial liquid velocity (U C ff = ds/dt - Uls), \ff is 
the effective area of the bubble on which the drag 
force acts; Sd is +1 or -1 for U e ff < 0 and U C ff > 0 re- 
spectively and C M c is the added mass coefficient. Ex- 
pressions for C M c, C DW i A eff are given in our previous 
work 16 . 

Buoyancy and gas momentum flux always promote 
bubble detachment, while surface tension tries to pre- 
vent it. The inertia force has two components. Bubble 
inertia represented by the first term of the inertia force 
expression is always an attaching force. On the other 
hand, liquid inertia, represented by the second term, as 
well as the liquid drag can be either attaching or de- 
taching, depending on whether U e ff > 0 or U e rr < 0. 

Bubble generation occurs in two stages. During the 
first stage, defined as the expansion stage, the bubble 
grows radially due to incoming gas flux, however the 
bubble base remains attached to the nozzle. Bubble 
volume and growth rate during the expansion stage arc 
respectively written as Vb = (tc/6) Db (t) and ds/dt = Vi 
dD B (t)/dt. A balance of the attaching and detaching 
forces marks the end of expansion stage (D B = D Bc ) 
and the beginning of the second stage, namely the de- 
tachment stage. During the detachment stage, the bub- 
ble continues to grow in size, the bubble volume being 
V B =(7i/6) D Bc 3 + Q d t. The bubble moves away from 
the nozzle, but still remains attached to it via a neck 
region. The bubble center located at a distance Y from 
the nozzle tip, moves at a velocity d s/dt = d Y/d t. 
Bubble motion during the two stages is described by a 
balance of forces acting on it: 

F B + F 0 4- F m + Fd + F| = 0 
At the end of the detachment stage, the bubble departs 
the nozzle base due to neck pinch-off. The neck 
pinches off when its length becomes equal to the noz- 
zle diameter. Therefore, the detachment criterion is 
written as: 

L N = Y -V4Db >D n 

Solving the force balance equations at the two stages 
of bubble formation, subject to the detachment crite- 
rion leads to a non-dimensional functional expression 
for the detached bubble diameter: 

D b * = f (Re P , We P , Fr P , U GS *, D N *,p*) 

Various dimensionless parameters are Reynolds num- 
ber Re P = P c UlsDp/|I c ; Weber number We P = 
p c U L s 2 D P /a; Froude number Fr P = p c Uls 2 /(p c - Pd)§ F> P ; 
Dimensionless superficial gas velocity ^ U G s = 
U gs /Uls; Dimensionless nozzle diameter D N = D N /D P 
and density ratio p* = p d /p c - Further details on the 
method of solution and predictions resulting from this 
theoretical model are discussed by Bhunia et. al. . 


A comparison of the present reduced gravity ex- 
perimental data with predictions of the numerical 
model is shown in figure 11. This figure displays 
variation of dimensionless bubble diameter with re- 
spect to non-dimensional superficial gas velocity. Two 
different sets of dimensionless nozzle diameter (D N ) 
and for each D N * two different Reynolds number con- 
ditions are considered. The numerical predictions show 
good agreement with the experimental data. Over a 
wide range of Re P , the present computational model 
predicts bubble diameter within ±10% of the experi- 
mental results. 
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Fig. 11: Comparison of numerical predictions with 
experimental results. Solid and dotted line - Numerical 
predictions, symbols - experimental data. 
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Fig 12: Comparison of computed dimensionless 
bubble formation time with reduced gravity experi- 
mental data. Solid and dotted line - numerical predic- 
tions, Symbols - experimental data; D N =0.1, Re P = 
2667, Wc P = 3.8, Fr P = 4.4; D N ’ = 0.2, Re P = 2318, 
We P = 3.9, Fr P = 8.0. 

The theoretical model is further validated by a 
comparison of the dimensionless bubble formation 
time (t* = t Uls/Dp) obtained from the model and the 
experimental data shown in figure 12. At low values of 
U G s\ the bubble formation time decreases sharply with 
increasing U G s , until it reaches an asymptotic limit of 
constant time, irrespective of the U GS value. It is fur- 
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thcr observed that with increasing nozzle diameter 
when bubble size increases, it takes longer for the 
bubble to detach. 

At higher superficial liquid velocity, the bubble de- 
viates from its spherical shape, as assumed in this 
model. The present model agrees well with the ex- 
perimental data up to a maximum Weber number of 
30, which corresponds to a Reynolds number of 7500 
for an air-water system using a 2.54 cm diameter pipe. 

CONCLUSIONS 

The present work focuses on bubble generation via 
single nozzle injection in a co- and cross-flowing liq- 
uid. This study is based on empirical data, obtained by 
performing experiments aboard the DC-9 Reduced 
Gravity Research Aircraft in parabolic flight. Effect of 
the flow conditions and geometry on detached bubble 
diameter and thereby void fraction of the resulting 
two-phase flow is investigated. It is shown that bubble 
diameter and void fraction increase with volumetric 
gas flow rate, pipe diameter and nozzle diameter, while 
they decrease with surrounding liquid flow. The im- 
portant role of the continuous liquid flow in detaching 
bubbles under reduced gravity conditions is thus em- 
phasized. 

It is of interest to note that bubble size and corre- 
sponding void fraction are somewhat smaller for the 
co-flow system than for the cross-flow configuration at 
similar flow conditions and flow geometry. From em- 
pirical evidence it is shown that the void fraction can 
be readily controlled in case of single nozzle gas in- 
jection by varying the flow geometry or the flow con- 
ditions. 

A theoretical model based on an overall force bal- 
ance acting on the bubble during the two stages of 
generation is also developed. Two sets of forces, one 
aiding and other inhibiting bubble detachment are 
identified. The theoretical model predicts bubble di- 
ameter in good agreement with the reduced gravity 
experiments. 
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